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Abstract 

We obtain exact moving and stationary, spatially periodic and localized solutions of a generalized discrete nonlinear 
Schrodinger equation. More specifically, we find two different moving periodic wave solutions and a localized moving 
pulse solution. We also address the problem of finding exact stationary solutions and, for a particular case of the 
model when stationary solutions can be expressed through the Jacobi elliptic functions, we present a two-point map 
from which all possible stationary solutions can be found. Numerically we demonstrate the generic stability of the 
stationary pulse solutions and also the robustness of moving pulses in long-term dynamics. 



I. INTRODUCTION 



tously [lj] 



The discrete nonlinear Schrodinger (DNLS) equation occurs ubiquitously [1[ throughout modern science. Most 
notable is the role it plays in understanding the propagation of electromagnetic waves in glass fibres and other optical 
waveguides as well as in the temporal evolution of Bose-Einstein condensates One of the variants of the DNLS 
model is the celebrated Ablowitz-Ladik (AL) model [J] which is an integrable model. Another aspect which stands 
out in favor of the AL model is that, while most other discrete DNLS models have stationary wave solutions [5[, 
this model has moving wave solutions. Further, these moving waves avoid the discreteness energy barrier (the so 
called Peierls-Nabarro (PN) barrier). These solutions have played a major role in the computational studies of 



the corresponding continuum NLS model 



(]] as well as in developing perturbation techniques Q] 



as well as in developing perturbation techniques [7j . It is clearly of 

10] and to try to obtain exact moving wave 



great interest to consider different variants of the DNLS equation 
solutions . The existence of such solutions might help in discovering new integrable models and would also 

help in further developing perturbative techniques in DNLS- type equations. The purpose of this paper is to report 
on the existence of exact moving as well as stationary solutions in a generalized DNLS model with seven parameters. 
For finite lattices we find two different periodic moving wave solutions while for the infinite lattice we find a localized 
moving pulse solution. 

In a recent paper, Pelinovsky Q has addressed the question of spatial discretization of the NLS equation with 
cubic nonlinearity 



ill + u xx + 2|u| 2 u = . (1) 
While the standard choice for the DNLS equation is 

iu n + u n+ i + it n _i - 2u n + 2\u n \ 2 u n = , (2) 

strictly speaking, there is no unique choice. Perhaps the only constraint on the corresponding discrete model is that 

in the continuum limit it should go over to the NLS Eq. (TTl) . By demanding that the semi-discretization is symplectic 

ft 

and few other requirements, Pelinovsky 9] showed that if one writes the DNLS equation in the form 

iu n + u n+ i + u„_i - 2u n + f(u n -i,u n , Un+i) = , (3) 

then the most general form for the nonlinear function / is given by 

/ = ai\u n \ 2 u n + a 2 \u n \ 2 (u n+ i + u n -i) + a 3 ul (w„+i + u n -i) 

+a4U n (\u n+1 \ 2 + |u„-i| 2 ) + a 5 u n (u n+ iu n ^i + u n -iu n+1 ) 

+a 6 U n (u^ +1 + + «7UnWn+lW«-l + a 8 ( I u n+ i 1 2 u n+ i + |lt„_i| 2 U n _i) 

+a g (u n -iu 2 l+1 + u„+it4_i) + aio(|u„+i| 2 u„_i + |u„_i| 2 u„ + i) , (4) 
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where u represents complex conjugate and the real valued parameters (et\, «io) satisfy the continuity constraint 

«i + a 7 + 2(a 2 + «3 + ck4 + a 5 + a 6 + a 8 + «9 + ctio) = 2 . (5) 

The purpose of this paper is to obtain moving as well as stationary solutions in this generalized model and study 
their stability. 

We note in passing that, under weaker constraints than that used in [9(, one can add to ([4]) the term proportional 
to u n (\u n -iu n \ + |u„u n+ i|), which was demonstrated to be translationally invariant and conserving the norm, E|u„| 2 

y 

The paper is organized as follows. In Sec. [Til we derive exact moving solutions for a seven-parameter DNLS model 
of Eq. ([¥]) with a± — ag — under the constraint ([5]). In addition, for a five-parameter translationally invariant 
DNLS equation we obtain a nonlinear map from which all possible stationary solutions can be derived. In Sec. IIIII 
we present numerical results for the stationary and moving pulse solutions to demonstrate their stability. Section HVl 
summarizes our main findings and concludes the paper. In the Appendix we list the identities for the Jacobi elliptic 
functions used in the derivation of the periodic wave solutions. 

II. ANALYTICAL RESULTS 

We now show that two moving periodic wave solutions can be obtained with this general cubic polynomial in case 
terms of the type |u„| 2 w„ are absent, i.e. 

oi\ = cts = • (6) 

It may be added here that the famous AL moving wave solutions are obtained in case only a-2 is non-zero while all 
other ai are zero. 

A. Solution I 

In particular, it is not difficult to show that one of the exact periodic wave solution to Eq. © [with / being given 
by Eq. (UJ) satisfying constraints © and ©] is given by 

u n = v4exp[— i(uot — kn + 6)] dsx[/3(n — v t + c),m] , (7) 

provided the following six relations are satisfied 



vf3 = 2A 2 (a 2 — a 3 ) sin(fe)cs(/3, m) , 



(8) 
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cs(/3, m)a% sin(2fc) + [otg sin (3 A;) — aio sm(&)]cs(2/3, m) = , (9) 
s * n j-^ = (a 2 — as) sin(fc)cs 2 (/3, to) — a6 sin(2fc)ds(/3, m)ns(/3, m) 

-[a 9 sin(3fc) - aio sin(fc)][cs 2 (2/3, m) + ds(2/3, m)ns(2/3, m)] , (10) 

[a 4 + a 6 cos(2fc)]cs(/3, to) + [ag cos(3fc) + aio cos(fc)]cs(2/3, m) = , (11) 
cos ( k } 

— j2 — = [ a 2 + a 3 ] cos(fc)cs 2 (/3, to) — [ai + a§ cos(2fc)]ds(/3, TO.)ns(/3, to) 
+ [2a 5 cos(2fc) + a^]cs((3, to)cs(2/3, to) 

— [ag cos(3fc) + aio cos(fc)][ds(2/3, to)iis(2/3, to) — cs 2 (2/3, to)] , (12) 
■^2 - -Jj = -2[a 2 + a 3 ] cos(fc)ds(/3, m)ns(/3, m) 

+2[a 4 + a 6 cos(2fc)]cs 2 (/?, m) - [2a 5 cos(2/c) + a 7 ]cs 2 (/3, to) . (13) 

Here c and <5 are arbitrary constants, k, lo, and v denote the wavenumber, frequency and velocity, re- 
spectively, of the periodic wave whereas cs(a, to), ds(a, to), ns(a, to) stand for the Jacobi elliptic functions 
cn/a, 77i)/sn(a, to), dn(a, TO.)/sn(a, to), l/sn(a, to) respectively with to being the modulus parameter (0 < to < 1) 



12j . While deriving these relations, use has been made of the local identities ([53")) to (|S"9")) for Jacobi elliptic functions 
dn(x, to) [13j which are given in the Appendix. 

It may be noted that Eqs. (|8|) to (|13j) determine the five parameters A, co, k, v, (3 and give us one constraint between 
the eight parameters a 2 , ■■■,aio (except as). In view of the constraint ([5]) between these parameters, it then follows 
that we have obtained a moving periodic wave solution with six parameters. As expected, in the limit a 2 7^ while 
all other a^ = 0, we recover the well known periodic wave solution of the AL problem Notice that in order that 
the periodic solution be compatible with the lattice, the modulus to has to be chosen such that (3N p = 2K(m) where 
K(m) denotes the complete elliptic integral of the first kind [2] and N p is the periodicity of the lattic 



B. Solution II 

As in the AL case, there is another periodic wave solution to the DNLS Eq. ^ with / being given by Eq. (|U) 
satisfying constraints ([5]) and ([6]). It is given by 

u n = A^/me~xjp[—i(Lut — kn + 5)] cn[/3(n — vt + c), to] , (14) 

provided the following relations are satisfied 

V0 = 2A 2 {a 2 - a 3 ) sin(fc)ds(/3, to) , (15) 
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«6 sin(2fc)ds(/3, to) + [ag sin(3fc) — aio sin(A;)]ds(2/3, m) = , (16) 

■rg ^ = (a 2 — CK3) sin(fc)ds 2 (/3, m) — a@ sin(2fc)cs(/3, m)ns(/3, m) 
-[a 9 sin(3fc) - a w sin(fc)][ds 2 (2/?, to) + cs(2/3, m)ns(2/3, m)] , (17) 

[d£4 + «6 cos(2/c)]ds(/3, to) + [ag cos(3fc) + aio cos(fc)]ds(2/3, to) = , (18) 

COS ( rZ) 

— P5 — = [ot2 + 013] cos(fc)ds 2 (3, to) — [a 4 + a§ cos(2fc)]cs(/3, m)ns(/3, to) 
+ [2a 5 cos(2fc) + a 7 ]ds(/3, m)ds(2/3, to) 

— [ag cos(3fc) + aio cos(fc)][cs(2/3, m)ns(2/3, m) — ds 2 (2/3, to)] , (19) 

- = -2[a 2 + a 3 ] cos(fc)cs(/3, m)ns(/3, m) 
+2[a 4 + a 6 cos(2/c)]ds 2 (/?, m) - [2a 5 cos(2/c) + a 7 ]ds 2 (/3, m) . (20) 

While deriving these relations, use has been made of the local identities (fuT)|) to (f6"6"|) for the Jacobi elliptic function 
cn(ir, to) [13! which have been given in the Appendix. 

As with the first solution, we again have a moving periodic wave solution with six parameters and again in the 
limit when only a 2 7^ while all other a.; are zero, we recover the well known periodic wave solution of the AL 
problem [l4|. In addition, note that in order that the periodic solution be compatible with the lattice, the modulus 
to has to be chosen such that f3N p — AK{m) where N p is the periodicity of the lattice 

C. Two- Point maps for stationary solutions 

With the ansatz u n (t) — /„e _IIJ ' we obtain from the DNLS Eqs. ©, (|U) the following second-order difference 
equation for the amplitudes 

/„_! - (2 - uj)f n + /„ +1 + a x fl + (a 2 + a 3 )/ 2 (/„-i + f n+1 ) + (a 4 + a 6 )/„ (f^_ t + / 2 +1 ) 

+ (2a 5 + a 7 )/«-i/«/n+i + a 8 {fLi + fn+l) + ("9 + Qio)/„-i/n+l(/n-i + fn+l) = 0. (21) 
For the following choice of parameters [that already includes the continuity constraint |5])] 

ol\ = a§ = 0, a 4 = — ag, &9 = — c^io, and a 7 + 2[a 2 + a.3 + 05] = 2 , (22) 
we get from (|2~Tj) the following second-order difference equation for the amplitudes 

/ n _l - (2 - w) /„ + /„ +1 + (a 2 + a 3 ) / 2 (/„_! + f n+1 ) + (2a 5 + a 7 ) /„-i/„/„+i - 0. (23) 



In this case, the stationary problem is exactly solvable. Indeed, one can obtain the first integral of (f2U)) and present 
it in the form of a two-point nonlinear map 

Zf n ± y/R(fn) 



fn+1 = (2-U))- 



Y 



2-ui + Yfl ' 
.(K-Xfi + ft), 



where 



R(fn) 



(2 - iof - K (2a 5 + a 7 ) 
2K (a 2 + a 3 ) (2a 5 + a 7 ) + 2 (2 - w) ' 
y = 2 (a 2 + a 3 ) Z + (2a 5 + a 7 ) , 

/fy 2 + (2-w) 2 (i-z 2 ) 

A. — 



(24) 



(25) 



(2-w)y 

Apart from the model parameters oti and frequency lo, the nonlinear map (|24p . (|25p contains the integration constant 
if. Due to the symmetry of equation (|2"5)) one can substitute for f n -i in (j24[) . For any set of admissible values 
fo, K, and w one can find the amplitudes of a stationary solution by iterating (f24l) . For R{f n ) > the map (|24|) 
gives two values for f n +i and one should take the one which satisfies the original three-point problem (|2"5|) . It is 
sufficient to take f n +i different from f n -i- 

The above two-point map can also be constructed from the Jacobi elliptic function solutions ||7J or (jT4j) as described 
in our recent work on a discrete 4 model \+\ . The corresponding DNLS equation has five free parameters because 
(|22p sets up five constraints between the ten parameters (o^) of the model. We note that any stationary solution to 
the DNLS equation defined by (|3|) and (j4|) with the parameters satisfying (|22p can be constructed from the nonlinear 



map 

equation 



(1251) . Such investigations have been carried out in our recent work on the DNLS equation 
15 



llj and the 



It is also worth pointing out that the three-point problem given by Eq. (|23p and the three-point problem studied 
by Quispel et al. [17| both can be presented in the following general form 

hl(fn) ~ h 2 (fn)fn-l 



f, 



n+1 



(26) 



M/n) _ M/n)/n-l 

For a particular choice of the functions hi(f n ), Quispel et al. have found a two-point map (i.e., the first integral of the 



corresponding three-point problem) which is quadratic in both <p n and n +i 



17j . For our choice of these functions, 



hi(fn) = (2 - w)/„ , h 2 (f n ) = 1 + ("2 + a 3 )f% , h 3 (f n ) = -(2a 5 + a 7 )f n , 



(27) 



we found the map (|24[) which is, in general, quartic in <p n and quadratic in 4> n +i- Clearly, our map (|24p does not 



belong to the 12 parameter map discussed in 



17|. 



The above result is new in that it generalizes the map reported in our recent work [151 ]. For completeness, let us 
also reproduce here the well-known result flQ for the case of 

a s — ag + aio, ol\ = «4 + chq, a.\ = 2a$ + a-j , and 4ai + 2[a 2 + a 3 + 2a$] = 2 , (28) 

when the continuity constraint ([5]) is satisfied and (f2Tj) reduces to the following second-order difference equation 



/ n _! - (2 - W )/„ + + a X /„ + + + J„-l/„+l] 

+ (a 2 + a 3 )fXfn-l + fn+l) + a S + + /„_!/„+!(/„_! + / n+1 )] = 0. (29) 



The first integral of |J29]) is 



V (/„_!,/») = + - (2 - U) f n -lf n + a, (fl_, + fl) /„_!/„ 

+ (a a + a 3 ) /n-x/" + "8 (fti + ft) +K = 0, (30) 



where K is the integration constant. This is so because (|29[) can be rewritten in the form 

V(fn, fn+l) — V(fn-1, fn) _ q 
fn+1 fn—1 

and one can verify that if V (/n-i, /n) = then (|29|) is satisfied. Solving the algebraic problem (|30|) iteratively for an 
admissible initial value /o one can construct a stationary solution to (|29p . This model has six free parameters because 
(|28p sets up four constraints between the ten parameters (a^) of the model. In general, stationary solutions to the 
DNLS equation with the parameters satisfying (|28[) cannot be expressed in terms of the Jacobi elliptic functions, 
but, as it was already mentioned, they can be constructed iteratively from (|30|) and they can be placed anywhere 
with respect to the lattice sites. 

We note that the translationally invariant discrete models possessing the form of equation (I31| have been introduced 
by Kevrekidis in 

D. Moving and stationary pulse solution 

In the limit m — > 1, both the periodic moving wave solutions {7} and (|14p reduce to the localized moving pulse 
solution 

u n = Aexp[— i(u>t — kn + 6)] scch[/3(rt — vt + c)] , (32) 

and the relations ([S]) to (fTB|) [as well as (fTS)) to (|2ll|) ] take a simpler form 

2sin(fc)sinh(/3) 

v= - p , (33) 
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2«6 sin(2fc) cosh(/3) + ag sin(3fc) — aio sin(fe) = , (34) 

[sinh 2 (/?) + (a 3 - a 2 ) A 2 ] sin(fc) = , (35) 

A 2 = 2 sinh 2 (/3) coshQ3) cos(fc) 

2 (a2 + CK3) cos(fc) cosh(/3) + 2 (ag — ae) cos(2fc) + a 7 — 2a 4 ' 

2[a 4 + a 6 cos(2fc)] cosh(/3) + ag cos(3fc) + aio cos(fc) = , (37) 

u = 2[1 - cos(fc) cosh(/3)] . (38) 

From (|33p . pulse velocity is zero when k — or fc = n. In the former case we have the non-staggered stationary 
pulse while in the latter case we have the staggered pulse. Needless to say that these remarks equally apply to the 
periodic wave solutions ((7|) and Q14p . In particular for k — we obtain the non-staggered stationary pulse solution 

u n = Aexp[— i(ut + 5)] sech[/3(n + c)] , 



A 1 



uj = 2 - 2cosh(/3), 
2 sinh 2 (/3) cosh(/3) 
2 (a2 + 03) cosh(/3) + 2 (a 5 — a 4 - a 6 ) + a.7 ' 

2 (a 4 + a 6 ) cosh(/3) + ag + aio = 0, (39) 



while in the latter case we obtain the staggered stationary pulse, 



A = 



Un = (~l) n Aexp[-i(ujt + 5)} sech[/3(n + c)] , 

w = 2 + 2cosh(/3), 
-2 sinh 2 (/3) cosh(/3) 
-2 (a 2 + a 3 ) cosh(/3) + 2 (a 5 — a 4 — a 6 ) + a 7 ' 

2 (a 4 -I- a 6 ) cosh(/3) — a g — aio = 0. (40) 



III. ANALYSIS OF THE PULSE SOLUTION AND NUMERICAL RESULTS 

For given model parameters a,, the moving pulse solution (|32"|) - (j38|) is characterized by two parameters, (3 > 
and —it < k < n. As it can be seen from Q33[) and psp . the pulse velocity and frequency do not depend on model 
parameters while the pulse amplitude does, see (I36[) . Using Q33|) one can express uj in (|38p as function of v and 
/3. Also using (|38|) one can express the group velocity dcu/dk. The pulse solution exists for given (3 and fe if the 
right-hand side of is positive and if (|3"4"1) . (133)) , and (|3T)l can be satisfied together with the continuity constraint 
, where we assume ((6]). 
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FIG. 1: (a) Velocity v and (b) frequency u> of the pulse as functions of the wavenumber parameter k at fixed value of the 
other parameter, inverse width of the pulse, f3 = 1/2 (dashed lines) and [3 — 1 (solid lines). These functions are defined by 
(|33|) and (|38|) and they do not depend on the model parameters a;. Pulse velocity is zero at k = and k = n, the former case 
corresponds to the non-staggered stationary pulse (|39|) while the latter case to the staggered stationary pulse J40J. 

As for the stationary pulse solution (|39p or (|40[) . for given model parameters oti, the moving pulse solution is 
characterized by a single parameter (5 > 0. In general, as far as the model parameters are fixed, the parameter (3 of 
the stationary pulse is also fixed through the last equation in f|39[) or (|40[) . However, for = — ae and ag = — aio, this 
constraint disappears and (3 can change continuously within a domain where A 2 > 0. Recall that in this particular 
case the stationary pulse solution can also be constructed from the two-point map presented in Sec. Ill C[ for which 
one should set the integration constant K = 0. 



A. Different moving solutions 

Coming back to the moving pulse solution (|32p ~(|38 p . several comments are in order. 

1. The relations (f33|) and (|38[) are exactly the same as in the AL case It is indeed remarkable that the velocity 
v and the frequency lo in our case are identical to those in the AL model even though our model has eight 
nonlinear terms (with coefficients a,i to aio with as = 0) while AL has only one term with «2 = 1- It is 
amusing to note that these two relations have also been obtained by Pelinovsky and Rothos from an entirely 
different approach [l^, namely from the linear dispersion relation for the corresponding differential advance- 
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FIG. 2: Nonzero model parameters 0:3 and a-j and the pulse amplitude A as functions of the parameter k at fixed value of the 
other parameter /3 = 1/2 (dashed lines) and /3 = 1 (solid lines). These functions are defined by (|43|l . For /3 = 1/2 the solution 
exists (i.e. A is real) for jfc| < 1.11 while for f3 = 1 it exists for |fc| < 1.23. The velocity v and frequency uj of the pulse are 
shown in Fig. [1] 




FIG. 3: Pulse amplitude A as a function of the parameter k at P = 1/2 (dashed lines) and /3 = 1 (solid lines) for the model 
with three nonzero parameters, 03, as and otj. Here we set 05 = 1 and find other model and pulse parameters from (144)1 . For 
P = 1/2 the solution exists (i.e. A is real) for |fc| < 7r/4 and 1.12 < |fc| < 37r/4, while for (5 = 1 it exists for |fc| < 7r/4 and 
1.24 < |fc| < 37r/4. The velocity v and frequency u) of the pulse are shown in Fig. [T] 
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FIG. 4: Pulse amplitude A as a function of the parameter k at /3 = 1/2 (dashed lines) and /3 = 1 (solid lines) for the model 
with three nonzero parameters, ct2, 0:3 and 05. The relation between model and pulse parameters are given by (|45[) . For 
/3 = 1/2 the solution exists (i.e. A is real) for |fc| < 1.26 and n/2 < \k\ < n, while for = 1 it exists for |fc| < 1.32 and 
7r/2 < |fc| < 7r. The velocity v and frequency u> of the pulse are shown in Fig. [T] 

delay equation. In Fig. [1] we show how v and o> depend on one of the pulse parameter, k, at fixed values of the 
other parameter, (3 = 1/2 (dashed lines) and (3 — 1 (solid lines). 

2. Unfortunately, we do not know the Hamiltonian from which the DNLS Eq. |3]) with / given by Eq. (j4]) can 
be derived. As a result, we cannot demonstrate the absence of the PN barrier from the energy consideration. 
However, since our stationary solutions have an effective translational invariancc (i.e. the solution is valid for 
any value of the constant c), this suggests that the PN barrier would be zero for these solutions. 

3. From Eq. (|35|) it follows that the moving pulse solution exists only if ct2 and/or 123 are nonzero. Further, in 
case «2 = 0, then it follows from Eq. ([33)) that < 0. 

4. In the limit when only ct2 is nonzero while all other on are zero, we recover the well known AL moving pulse 



solution 



4]- 



5. The sn-type and hence dark soliton solution can also be obtained in this generalized model provided the right 
hand side of the continuity Eq. (0) is —2 (instead of 2). 

6. In case only ai and/or 013 are nonzero while all other = 0, then the generalized DNLS Eq. ([3]) with / given 
by Eq. ((¥]) conserves the momentum defined by 

P = i ^ { U n+lUn ~ U n+1 U n ) . (41) 

n 

On the other hand, in case only a§ and/or a*j are nonzero while all other ai = 0, then the generalized DNLS 



12 

Eq. §3§ with / given by Eq. ^ conserves the momentum defined by 

P = i ^ ( U n+2U n - u n+2 u n ) . (42) 

n 

Expression (|4*2"j) is similar to that introduced in [15| for the c/> 4 discrete equation. 

7. From Eqs. to it follows that the moving pulse solution is also possible when only two of the eight 

parameters are nonzero. For example, the moving pulse solution (|32p exists in case 03,07 are nonzero while 
all other cti are zero. While the relations ((33)) and (|38|) are always valid, the other relations and the constraint 
§5§ take the form 

= 1 =2f1- ^ A 2 - C0 K fc ) C0S Mff)sinh 2 (/3) U on 

3 l-2cos(fc)cosh(/3)' 7 1 3j ' o 3 [cos(fc)cosh(/3)-l] + l' [ 1 



For a pair of pulse parameters, k and /3, we find 03 and then 07 and A from (|43|) and present the result in 
Fig. [2 for /3 = 1/2 (dashed lines) and /? = 1 (solid lines). For 0=1/2 the solution exists (i.e. A is real) for 
\k\ < 1.11 while for (3 = 1 it exists for |fc| < 1.23. One can see that the non-staggered stationary pulse (k = 0) 
exists while staggered stationary pulse (k = it) does not exist in this case. 

8. The moving pulse solution (|32[) also exists in case only (i) 013,05 are nonzero; (ii) 02,05 are nonzero and 
cos(2fc) = 0, i.e., regardless of the model parameters, in this case one can have only k — ±ir/4 and k — ±37r/4. 
Constraints similar to those in (|43|) are easily written down from relations (O and (|33|) to (|38|) . We were 
unable to find other sets of model parameters supporting the pulse solution when there are only two nonzero 
parameters. 

9. There are several possibilities, with three of the eight o^ being nonzero (the remaining five ai being zero), in 
which case the moving pulse solution (|32|> is still valid. These cases are: (i) 02, 03, 05 are nonzero; (ii) 02, 03, 07 
are nonzero; (iii) 03,05,07 are nonzero; (iv) 02,05,07 are nonzero; (v) 02,04,06 are nonzero with 04 = og 
and k = ±7r/2; (vi) 03,04,06 are nonzero with 04 = 06 and k — ±ir/2; (vii) 02,09,010 are nonzero with 
09 = oio and k = ±7r/4 or k = ±37r/4. 

In all these cases the constraints similar to those in (|4"3"]l are easily obtained from relations ([5]) and to (|38[) . 
For example, in case only 03,05,07 are nonzero, while the relations (133|) and (138p are always valid, the other 
relations and the constraint (J5]) take the form 



1 — 205 sin 2 (k) . ) A 2 cos (fc) sinh 2 (/3) cosh (/?) 

l-2cos(fc) cosh(/3)' 011 ~ { - a s-a 5 ), - - - — — ^ — — - - - ^ ^ 2 ^ . 



The number of constraints in this case is such that one has a free model parameter, say 05, and pulse parameters 
k and [3 can change continuously within a certain domain. For 05 with a small absolute value the solution is 
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close to (U2J) shown in Fig. [5J but, for example, for as = 1 the result is qualitatively different, as it can be seen 
from Fig. [3J Also note that in this case the non-staggered stationary pulse (k = 0) exists while the staggered 
stationary pulse [k — tt) does not exist. 

On the other hand, in case only «2, 0^3, a§ are nonzero we have the following constraints 

— a5Cos(2fc) 2 cos(fc) sinh 2 (/3) cosh(/3) 

2 cos(fc) cosh(/3) ' " " ' (a.2 + 0:3) cos(fc) cosh(/3) + a$ cos(2fc) 

The relation between pulse parameters and model parameters in this case are shown in Fig. 3J In this case one 
has both non-staggered and staggered stationary pulse solutions for k = and k — tt, respectively. 

We give two more solutions, for the case when only oli, ct4, a§ are nonzero, 

2 sinh 2 (3 1 — a2 tt 

A = , 04 = — - — , a 4 = a 6 , k = ±-, (46) 

ct2 2 2 

and for the case when only a 2 , ag, aio are nonzero, 

2 sinh 2 ^ 1 - a 2 , , 7r , , 37r 

A = , aio = — ^ — , a 9 = aio, fe = ±— , or fc = ± — . (47) 

tt2 2 4 4 

These two moving solutions are interesting because for them the relations (f22|) are violated. These models have 
one free parameter, for example, a 2 > 0. Among the two pulse parameters, only (3 can change continuously, 
while k can assume only a few isolated discrete values, that do not depend on model parameters en. For the 
cases when there are only three nonzero parameters, we were unable to find sets of model parameters supporting 
the pulse solution other than the ones described above. 

10. Similarly, there are several possibilities when less than eight parameters are nonzero and still the moving pulse 
solution (f3"2"]) continues to exist and relations similar to those in Eq. (|4"3"]l can easily be obtained in all these 
cases. 

11. Since the DNLS equation ([3]), ([4]) with any set of parameters ctj satisfying the continuity constraint ([5]) reduces 
to the same NLS equation |T]), for a sufficiently wide (small j3) and slow (small |fc|) pulse, all the solutions given 
above are close and can be well approximated by the moving solution to the continuous NLS equation. 



B. Stability of the pulse solution 



Let us now discuss the small amplitude vibration spectrum for the lattice containing a stationary pulse in order to 
observe the peculiarities of the spectrum of the pulse in a translationally invariant lattice and to discuss the stability 
of the pulse. The vibrational spectrum was calculated following the methodology presented in Ref. [20j similar to the 
work in [ll| . In brief, we consider a small complex perturbation of a stationary solution and substitute the ansatz 
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« n (t) = lf n + e n (t)]e~ luJt with e n (t) = a n (t) + ib n (t) into the DNLS equation ©, (0} and obtain a linear equation 
for e n (t). Separating real and imaginary parts of the equation, we derive the following system 

(48) 

where vectors a and b contain a n and b n , respectively, while the nonzero coefficients of matrices K and J are given 

by, 

K n , n -\ = 1 + («2 + a 3 ) fl + (2a 5 + a 7 ) f„f n +i + (ocg + a w ) (fn+i + 2/„-i/, t +i) , 
K n ,n = - (2 - lo) + 2 (a 2 + a 3 ) /„ (/„_! + f n+1 ) + (a 4 + a 6 ) (f^_ 1 + fl +l ) + {2a 5 + a 7 ) /„-i/„+i, 

K n!n+ i = 1 + (a 2 + a 3 ) fl + (2a 5 + a 7 ) /„_i/„ + (a g + a w ) (f„-i + 2/„-i/„+i) , (49) 

Jn,n-1 = -1 - («2 - Ct3)fn ~ 2a 6 /„-l/ ti - (X7 fnfn+1 ~ 2a 9 / n _i/ n+ i + (ag - a 10 ) fn+i, 
Jn,n = (2 - W) - 2a 3 /n (fn-1 + fn+l) - («4 " "e) (/n-1 + fl+l) - ( 2a 5 - OL 7 ) fn-lfn+1, 
Jn,n+1 = -1 - ("2 - a 3 ) - 2a 6 /n/n+l ~ "7/n-l/n ~ 2a 9 /n-l/n+l + («9 ~ «lo)/^-i- (50) 

A stationary solution is characterized as linearly stable if and only if the eigenvalue problem 

= 7(1 (51) 




results in nonpositive real parts of all eigenvalues 7. 

Setting in the above matrices f n — 0, and solving the resulting eigenvalue problem one finds the spectrum of 
vacuum 



= ± 



-u + 4 sin" 2 I — 



(52) 



where Vt and Q are the frequency and the wavenumber of a small-amplitude harmonic mode, respectively. A stationary 
pulse was placed in the middle of a lattice of N = 200 points and the eigenvalue problem (|5Tj) was solved employing 
periodic boundary conditions. Here we do not aim to present a comprehensive numerical study of the stability of 
the pulse because the DNLS equation under consideration has a multi-dimensional parameter space and such an 
exhaustive study would entail enormous effort. Instead, our intent is to check several sets of parameters and to 
provide a few examples illustrating the generic stability of the pulse solution. 

Two examples of stationary, stable pulses and their spectra are presented in Fig. [5j Left panels present the 
results for a non-staggered pule, while the right panels are for a staggered pulse. Model parameters correspond to 
a translationally invariant lattice, i.e., they satisfy (|2"2"]) . For panels (a), (b) parameters are a 2 = 1, a 3 = —1/2, 
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«4 = —as = 1/2, «5 = —1/2, «7 = 2, and ag = — aio = —1/2. For panels (a'), (b') parameters are a,i = 1, = 0, 
a 4 = — a e — 1/2, a 5 = —1/2, 017 = 1, and ag — — a 10 = —1/2. The choice of parameters is rather arbitrary. For 
the non-staggered pulse all coefficients are nonzero so that all terms of the DNLS equation are involved. For the 
staggered pulse we found that nonzero makes the pulse unstable, that is why we set this coefficient equal to zero. 
The non-staggered and staggered pulses are defined by, respectively, ((39)) and (|40[) with parameters /3 = 1, S = 0, 
and c = 0.25. We then found for the non-staggered pulse uj = —1.0862 and A = 1.2946, and for the staggered one, 
u> = 5.0862 and A = 1.1752. The pulses are placed asymmetrically with respect to the lattice, nevertheless, they 
are stationary and stable since all eigenvalues 7 have zero real parts. The spectrum of non-staggered pulse contains 
the spectrum of vacuum ([52")) with the bands 1.0862 < < 5.0862; the three pulse internal modes with frequencies 
±0.195, ±5.290, and ±6.396; and the two pairs of zero eigenvalues, one pair corresponding to the translational 
invariance and another to the invariance with respect to the phase shift. The spectrum of staggered pulse is similar 
but it contains not three but only one pulse internal mode with frequencies ±5.308. 

We have checked the stability of stationary pulses (both non-staggered and staggered) with different (3 ~ 1, and 
also different positions with respect to the lattice, c, and for various model parameters with \a.i\ ~ 1, and in many 
of the cases found these pulses to be stable. Thus, we conclude that the stationary pulse solutions (f3"9"| and (|4"0"|) to 
DNLS equation ©, (0| with parameters satisfying (f2"2"| are generically stable. 

We have also checked the stability of a stationary pulse in the model where the pulse solution exists only for a 
selected (3 and, for the pulse placed asymmetrically with respect to the lattice we found that it is stable. In this 
simulation, for the solution (|39|) we took the following pulse parameters f3 = 1/2, c = 1/4, and model parameters 
u-2 = —0.2553, 04 = a6 = —1/2, and ag — 2.2553 with all other on equal to zero. We then found ui = —0.2553 and 
A = 0.6557. 

The robustness of moving pulse solutions was checked by observing the evolution of their velocity in a long-term 
numerical run. For pulses with amplitudes A ~ 1 and velocities v ~ 0.1 and for various model parameters supporting 
the pulse, \ai\ ~ 1, we found that the pulse typically preserves its velocity with a high accuracy. Two examples of 
such simulations, one for the non-staggered pulse and another one for the staggered pulse are given in Fig. [6] (a), 
(b) and (a'), (b'), respectively. In (a) and (a') we show the pulse configuration at t = Q and in (b) and (b') the 
pulse velocity as a function of time for two different integration steps, r = 5 x 10~ 3 (solid lines) and r = 2.5 x 10~ 3 
(dashed lines), while the numerical scheme with an accuracy 0(t ) is employed. 

In both cases, one can notice the linear increase in the pulse velocity with time, which is due to the numerical error, 
since the slope of the line decreases with the decrease in r. The presence of perturbation in the form of rounding 
errors and integration scheme errors does not result in pulse instability within the numerical run. Velocity increase 
rate for the staggered pulse in (b') is larger than for the non-staggered one in (b). This can be easily understood 
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because the frequency of the staggered pulse in almost five times larger than that of the non-staggered one. 

The pulse presented in Fig. [5] (a) is given by (|44p. The model has one free parameter and we set ccg = 1. For 
the pulse parameters we set — 1 and k = 0.102102 (close to zero). Then we find from (|32|) . (|33|) . ((38|) . and (|44|) 
the pulse velocity v = 0.239563, frequency w = —1.07009, and amplitude A = 1.7087, and the dependent model 
parameters a 3 = -0.473034 and a 7 = 0.946068. 

In Fig. [6] (a') the moving pulse solution is given by (|45|) . The model has one free parameter and we set = 0.3. 
For the pulse parameters we set = 1 and k = 3.09447 (close to ir). Then we find from ([52^1. (j3"3"|) , and (|4"5]) 

the pulse velocity w — 0.110719, frequency u> — 5.08274, and amplitude A — 1.65172, and the dependent model 
parameters a 2 = 0.603116 and a 3 — 0.0968843. 

Similar results were observed for the cases when only a 3 and 0:5 are nonzero; only a 3 and a.7 are nonzero; only 
Q!2, 013, and «5 are nonzero; only 012, a 3 , and ay are nonzero; and only a 2 , ct5, and are nonzero. 

So far we have studied numerically the pulses in the models with the parameters satisfying (|22[) . However, moving 
pulse solutions exist even in the case when (f2"2"| is violated. Two such solutions are presented by (|4"6")l and (|4"T|) 
together with (f3"2")l . (|33|) . and (f3"5)) . As it can be seen from Fig. the pulses show a stable long-term dynamics with 
pulse velocity being practically constant with the accuracy increasing with decrease in the step size of numerical 
integration. The pulse in (a) is given by (f3"2")) , (f3"3")) , (f3"5| . and (|4"6"]) . The model and pulse parameters are as follows: 
a 2 = 2, a 4 = -1/2, and a 6 = -1/2; /3 = 1, k = tt/2, w = 1.662, w = -0.1822, and A = 0.8310. The pulse in (a') 
is given by (|32)) , (|33|) . ([38| . and (|47)) . The model and pulse parameters are as follows: a 2 — 2, ag = —1/2, and 
otic = -1/2; = 1, k = tt/4, u = 2.350, w = 2, and A = 0.8310. 

Velocity increase rate in (b) is considerably larger than in (b') (note the different abscissa scale for these two 
panels) and this result can be expected when we take into account that pulse frequency in (b) is 11 times larger than 
in (b'). 

IV. CONCLUSIONS AND FUTURE CHALLENGES 

For the nine-parameter DNLS equation ©, (U|) with the continuity constraint (O, in Sec. Ill Al and Sec. Ill B[ 
we obtained the two moving periodic wave solutions for the case of ot\ — as = (thus, the moving solutions are 
supported by the seven-parameter model). The solutions have the form of dn and cn Jacobi elliptic functions. In the 
limit m — > 1 both solutions reduce to the moving pulse solution (see Sec. Ill D[) . We found and described several sets 
of model parameters supporting the moving pulse solution. For the particular choice of model parameters (|22j) . the 
problem of finding stationary solutions is integrable and the first integral of this problem was given in Sec. Ill Cl in the 
form of a nonlinear map. From this map any stationary solution of the corresponding problem can be constructed 
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FIG. 5: Two examples of (a), (a') stationary pulse profiles and (b), (b') their spectra. Left panels show the results for a 
non-staggered pulse, while right panels are for a staggered pulse. Model parameters correspond to a translationally invariant 
lattice, i.e., they satisfy (|22p . For (a), (b) parameters are a?2 = 1, as = —1/2, Q4 = — = 1/2, a 5 = —1/2, 07 = 2, 
and ag = — Q10 = —1/2- For (a'), (b') parameters are Q2 = 1, 03 = 0, 04 = — as = 1/2, 0:5 = —1/2, 07 = 1, and 
ag — — aio = —1/2. The pulses are denned by, respectively, (139[) and (|40|l with parameters /3 = 1, 5 = 0, and c = 0.25. Pulses 
are placed asymmetrically with respect to the lattice, nevertheless, they are stationary and stable since all eigenvalues 7 have 
zero real parts. The spectra also contain two pairs of zero eigenvalues, one pair corresponds to the translational invariance 
and another to the invariance with respect to the phase shift. 



iteratively. 

We found the stationary pulse solutions to be generically stable, i.e., for rather arbitrary choice of model parameters 
|ct;| ~ 1, in many cases, the spectra of the small-amplitude vibrations calculated for the lattice containing a pulse 
included no eigenvalues with positive real parts. In addition, we confirmed the robustness of moving pulses by 
observing the pulse velocity evolution in a long-term numerical run. We found the velocity to be nearly constant 
and the deviation from constancy was attributed to the influence of the accuracy of the numerical integration. We 
specifically note that the moving pulse solutions exist and they exhibit a stable behavior in long-term numerical runs 
even for models which do not support translationally invariant stationary pulse solutions, as demonstrated in Fig. [7] 

On using the identities for the Jacobi elliptic functions cn and dn given below and similar identities for sn, one 
can similarly obtain exact solutions of a rather general discrete A</> 4 field theory with four parameters, as well as 



of a modified Fermi-Pasta-Ulam (FPU) model 



2l|, which will be discussed elsewhere. Our results are potentially 
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FIG. 6: (a) Non-staggered moving pulse at t = and (a') same for the staggered pulse. In (b) and (b') the long-term 
evolution of pulse velocity is shown for the corresponding pulses for the integration steps of r = 5 x 10 -3 (solid line) and 
t = 2.5 x 10~ 3 (dashed line). Numerical scheme with an accuracy 0(t 4 ) is employed. Pulses preserve their velocity with the 
accuracy increasing with the increase in the accuracy of numerical integration. Within the numerical run, the pulse dynamics 
is stable in spite of the presence of small perturbations in the system in the form of rounding errors and integration scheme 
errors. The pulse in (a) is given by (|32p . (|33p . (J38J) , and (|44[) . The model and pulse parameters are as follows: ay = —0.473034, 
a 5 = 1, and a 7 = 0.946068; f3 = 1, k = 0.102102 (close to 0), v = 0.239563, u = -1.07009, and A = 1.7087. The pulse in (a') 
is given by (|32|l . (|33p . (|38|l . and (|45p . The model and pulse parameters are as follows: «2 = 0.603116, C13 — 0.0968843, and 
a 5 = 0.3; /3 = l,k = 3.09447 (close to tt), v = 0.110719, u = 5.08274, and A = 1.65172. 

important for optical pulse propagation in glass fibers and optical waveguides [2] and time evolution of Bose-Einstein 
condensates [aj. 
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FIG. 7: Results similar to that shown in Fig. [6] but for models that are not translationally invariant, (a) and (a') show the 
moving pulse profiles at t = 0. In (b) and (b') the long-term evolution of pulse velocity is shown for the corresponding pulses. 
The integration steps are (b) r = 10~ 3 (solid line) and r = 5 x 10~ 4 (dashed line) and (b') r = 5 x 10~ 3 (solid line) and 
t — 2.5 x 1CP 3 (dashed line). The pulse in (a) is given by (|32|l . (|33|l . (I38|l . and (|46l) . The model and pulse parameters are as 
follows: a 2 = 2, a 4 = -1/2, and a 6 = -1/2; /3 = 1, k = tt/2, v = 1.662, u = -0.1822, and A = 0.8310. The pulse in (a') 
is given by (j32|> . (|33p . (|38|l . and (|47p . The model and pulse parameters are as follows: ct2 — 2, ctg = —1/2, and aio = —1/2; 
p = l,k = tt/4, « = 2.350, uj = 2, and A = 0.8310. 

V. APPENDIX 

We list here the various identities for the Jacobi elliptic functions dn(a;, m) and cn(a;, m) which have been used in 
obtaining the various solutions in this paper. 
Identities for dn(x,m) 

dn 2 (ir, m)[dn(x + a, m) + dn(x — a, to)] = — cs 2 (a, m)[dn(a: + a, to) +dn(x — a, to)] + 2ns(a, m)ds(a, m,)dn(x, to) , (53) 

dn(x, TO)dn(x + a, m)dn(x — a, to) = — cs(a, m)cs(2a, to) [dn(x + a, to) + dn(a; — a, to)] + cs 2 (a, m)dn(x, to) , (54) 

dn(x, TO)[dn 2 (ir + a, to) + dn 2 (x — a, to)] = ds(a, m)ns(a, m)[dn(x + a, to) + dn(x — a, to)] 

— 2cs 2 (a, TO,)dn(x, to) + tocs(o, m)[cn(x + a, m)sn(x + a, to) — cn(x — a, TO)sn(x — a, to)] , (55) 



dn(a; + a, m)dn(x — a, to) [dn{x + a, m) + dn(x — a, to)] 

= [ds(2a, m)ns(2a, to) — cs 2 (2a, to)] [dn(x + a, to) + dn(.x — a, m)] 

+TOCs(2a, to) [cn(a; + a, TO.)sn(:r + a, to) — cn(a; — a, TO)sn(x — o, to)] , 

dn 2 (x, TO)[dn(x + o, to) — dn(:r — a, to)] = — cs 2 (a, TO)[dn(x + a, to) — dn(:r — a, to)] — 2mcs(a, m)cn(x, m)sn(x, to) 

dn(x, TO)[dn 2 (x + a, to) — dn 2 (a; — a, to)] = ds(a, TO,)ns(a, TO)[dn(x + a, to) — dn(x — a, to)] 
+TOCs(a, to) [cn(x + a, to)sii(x + a, to) + cn(x — a, TO)sn(x — a, m)] , 

dn(x + a, TO,)dn(.x — a, to) [dn(.T + a, to) — dn(x — a, to)] 

= [ds(2a, TO)ns(2a, to) + cs 2 (2a, to)] [dn(x + a, to) — dn(x — a, to)] 

+TOCs(2a, to) [cn(£ + a, TO,)sn(:r + a, to) + cn(x — a, TO)sn(x — a, to)] . 

Identities for cn(x, to) 



men 2 (x, m) [cn(x + a, to) + cn(x — a, to)] = — ds 2 (a, to) [cn(x + a, to) + cn(x — a, to)] + 2ns(a, to)cs(o, m)cn(x, to) 

TOcn(x, TO)cn(a; + a, TO)cn(x — a, to) = — ds(a, m)ds(2a, TO)[cn(a; + a, to) + cn(x — a, to)] + ds 2 (a, m)cn(x, to) , 

mcn(x, TO)[cn 2 (x + a, to) + cn 2 (a; — a, to)] = cs(a, m)ns(a, m)[cn(x + a, to) + cn(x — a, to)] 
— 2ds 2 (a, m)cn(x, to) + ds(a, TO)[dn(x + a, TO)sn(a; + a, to) — dn(x — a, m)sn(a; — a, to)] , 

mcn(x + a, m)cn(x — a, to) [cn(x + a, to) + cn(x — a, to)] 

= [cs(2a, TO)ns(2a, to) — ds 2 (2a, TO,)][cn(a; + a, to) + cn(x — a, to)] 

+ds(2a, to) [dn(x + a, to)sii(:z! + a, to) — dn(x — a, TO,)sn(a; — a, to)] , 

tocii 2 (x, TO)[cn(x + a, to) — cn(x — a, to)] = — ds 2 (a, TO,)[cn(iz; + a, to) — cn(:r — a, to)] — 2ds(a, TO)dn(ir, m)sn(x, to) 

mcn(x, TO)[cn 2 (ir + a, to) — cn 2 (x — a, to)] = cs(a, m)ns(a, TO)[cn(a; + a, to) — cn(x — a, to)] 
+ds(a, TO)[dn(x + a, TO)sn(x + a, to) + dn(x — a, to)sii(.t — a, to)] , 
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mcn(x + a, m)cn(x — a, to) [cn(x + a, m) — cn{x — a, to)] 

= [cs(2a, TO)ns(2a, m) + ds 2 (2a, m)][cn(a; + a, to) — cn(x — a, to)] 

+ds(2a, to) [dn(x + a, TO)sn(a; + a, to) + dn(a; — a, TO)sn(a; — a, to)] . (66) 
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